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,-_i I Abstract 

_i I We discuss the stability of strangelets by considering dynamical chiral symme- 

^ ' try breaking. For quark droplets of finite volume, we formulate the Nambu-Jona- 

^«0 . Lasinio model with a basis set of the quark wave functions in the chiral bag model. 

r^ , Chiral symmetry breaking for the finite volume bag is discussed in a mean field 

Tij- I approximation. Effects of the pion cloud including the chiral Casimir effect are in- 

^^ ■ vestigated. Physical quantities of the quark droplets such as masses and radii are 

f— ^ . obtained for quark droplets of baryon numbers A < 5. We also apply our model 

"^ I setting to the nucleon and discuss the stability of the quark droplets against nuclei. 

Q^ 1 Introduction 

-f^. : 

>■ ! The study of the strange matter containing a large amount of strangeness is one of inter- 
'k> ! esting subjects in the hadron and quark physics. It has been considered that the strange 

^ I matter can exit to be a stable state due to a large number of degrees of freedom of the 
color and the strangeness [U El Ej- It was proposed that the strange matter could be 
the true ground state of the matter of the strong interaction. A finite volume droplet of 
the strange matter, which is called strangelet, has been also investigated. In the analysis 
of the MIT bag model, it was shown that the strangelets were stable particles of large 
baryon number and small electric charge as compared with normal nuclei jHl E] • 

Several scenarios of formation of the strangelets have been considered in various con- 
texts. The strangelets may be formed in the relativistic heavy ion collisions in the acceler- 
ator which can produce the quark-gluon plasma [310]. In the cosmic scale, the strangelets 

*yasui@th. phys.titech.ac.jp 
^hosakaQrcnp. osaka-u.ac.jp 



may be formed in the QCD phase transition in the early universe and/or in the colhsions 
of the compact stars containing a rich amount of strange quarks [7]. As a remnant of 
such processes in the universe, the strangelets may be observed as exotic particles in the 
cosmic rays. Indeed, there are several reports on the observation of the exotic particles 
in the cosmic rays [3 El El UHl UH 1121 HHj- The baryon number of the observed particle is 
from a few hundreds to thousands, while the electric charge is the order of ten, which are 
very much different from the ordinary nuclei. 

In the early stage of the study of the strangelets, it was assumed that the quark matter 
was described as almost free Fermi gas confined in a bag, in which chiral symmetry was 
restored [Slil]- There, the finite current mass of quarks played only a minor role, because 
the Fermi energy was larger than the current quark masses of the ud and s quarks. On the 
other hand, when chiral symmetry is maximally broken, the s quarks would not appear in 
the ground state of the quark matter due to the heavy mass in comparison with that of the 
ud quarks. By using an effective model of QCD, Buballa et al. showed that the strange 
matter was not the stable state of the quark matter at low density, while it could become 
the ground state at high density, more than a few times of the normal nuclear matter 
density ^M, IT^ IT^ . At high density, the chiral symmetry of the s quark is restored, and 
the s quark mass becomes smaller than the Fermi energy of the ud quarks. Therefore, 
the weak decay process is allowed for the generation of s quarks in the quark matter. 
However, this model analysis was not applied to the strangelets of finite volume, since the 
surface effects were not taken into account. 

For the discussion of the stability of the strangelets with dynamical quarks, the con- 
stituent quark model was used in Refs. [T71 UHl HHj • There, it was concluded that the 
strangelets were not stable objects due to the large s quark mass. In order to discuss 
the dynamical generation of the quark mass, the point-like interaction of the Nambu- 
Jona-Lasinio (N JL) type interaction was introduced in the quark droplet [201 1^ 1221 120] ■ 
There, the quark wave function of the MIT bag model was used as a set of basis functions 
in the finite volume system. In that discussion, it was shown that the chiral symmetry in 
the s quark sector was restored for the baryon number A < 10^. There, the strangelets 
were the ground state of the quark droplets. The masses of the strangelets were also 
obtained. It was concluded that the strangelets were not stable objects as compared with 
the normal nuclei, since the energy per baryon number E/A of the strangelets were found 
to be larger than the mass of the ground state baryons. 

In these previous studies, the pion and kaon cloud around the bag was not taken into 
account, though it was necessary to conserve chiral symmetry at the bag. For the quark 
droplets with relatively small baryon number A < 10, the MIT bag with a scalar type 
boundary condition |211 ESI I2ni I2Z| does not give a fully chirally symmetric formulations. 
In this paper, we discuss effects of the pion cloud on the chiral symmetry breaking in the 
quark droplets for small baryon number by using the chiral bag model. 

The contents of this paper are as follows. In Section 2, we introduce the NJL lagrangian 
with a basis set of the chiral bag model as an effective model of the quark droplet. In this 
paper we call this hybrid model as the NJL chiral bag model. In Section 3, we show the 
numerical results and discuss the stability of quark droplets. Section 4 is devoted to the 



conclusion. 



2 Model 



2.1 NJL chiral bag model 

We use the NJL type interaction as an effective theory for dynamical chiral symmetry 
breaking |2H]- For the construction of the quark wave function in a finite volume, we 
consider a bag model for the quark droplets. Here, we consider the chiral bag model to 
conserve chirality at the bag surface [23 EHl 1211 IHD] • Then the quarks inside the chiral 
bag are interacting through the four-point interaction of the NJL model, and the meson 
cloud exists outside the bag. The idea of a chiral bag model with the NJL interaction was 
first presented by T. Kunihiro in 1984 pTll32j^. In this paper, we call this hybrid model 
as the NJL chiral bag model. The model lagrangian is given by 
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where ip = {u, d, s)* is the quark field, and the current mass matrix mo = diag(moM, mod, rriQs 
The second term in the first bracket is the NJL type interaction term invariant under 
U{Nf)L X U{Nf)R with Nf = 3, in which A^ (a = 0, ■ ■ ■ , 8) is the Gell-Mann matrices 
normalized by trAaAf, = 26ab- We do not consider the 't Hooft term for explicit U{1)a 
breaking. The step functions are multiplied in order to distinguish the quark phase inside 
the bag and the meson phase outside the bag. Here, r is a distance from the center of 
the bag. We assume that the strangelet has a spherical shape with a bag radius R. The 
second term with the delta function realizes the chiral boundary condition for vector and 
axial vector current fiows at the bag surface with keeping chiral symmetry, in which we 
define 
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The third term is the meson lagrangian outside the bag, where the meson field is 
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The meson fields are normalized by the pion decay constant /^ = 0.093 GeV. We introduce 
the Skyrme term to obtain stable solutions [3311311, where the coefficient is set e = 6.17 
so as to reproduce -Eat+a = 1.1 GeV, the average mass of N and A. The meson mass 
term is given by 
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where M r^ m. In the following discussion, however, we turn off the pion mass term which 
does not change the result qualitatively. 

In the meson sector, we consider only the vr meson by neglecting K and r] mesons, 
because they are heavier than the vr meson. Therefore, the meson field is written as 
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where we define u.,^ = e^'^"" . 

In the quark sector, chiral symmetry breaking is caused by the non-perturbative ef- 
fect in the NJL type interaction, which gives a non-zero expectation value of the quark 
scalar condensate {qq) . In the NJL interaction, we consider the mean field approximation 
[qqY — > '^'i'ii'il) — {qq)'^ for q = u, d and s, respectively. We define the dynamical quark 
mass 
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for each flavor. Note that mQg = for q = u, d, but rrios 7^ for g = s. Then, the 
lagrangian in the quark sector in Eq. (^ is written as the sum of the ud and s quark 
sectors 



where Cud is the lagrangian in the ud quark sector 
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and Cs is the lagrangian in the s quark sector 
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In Eq. ©, we define u1^ = e^'^"^'^^\ The ud quarks interact with the pion of the hedgehog 
solution, while the s quarks do not. We define the hedgehog wave function h in the ud 
sector as discussed below. We set rrih = rriu = ma and rrioh = 0. 



2.2 Hedgehog solution 



We assume that the meson field and the interaction are strong, and therefore the non- 
hnear effects should be included to all orders. In order to consider the non-linear effect, 
we assume the hedgehog ansatz for the tt meson sector. The hedgehog ansatz conserves 
the grand spin K = J + 1 with total angular momentum J and isospin I. The hedgehog 
pion field for r > i? is then written as 
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with the chiral angle F{r) as a function of r and a unit radial vector n in the coordinate 
space. In the hedgehog ansatz, the unit radial vector is identified with that in the isospin 
space. Then, we obtain the energy of the pion field 
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Here, F in the left hand side is the chiral angle at the bag surface. The profile function 
F[r) is obtained by solving the equation of motion 
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which is given by taking a variation of the energy F^^ with respect to Fir) . The asymptotic 
form of this solution in the limit of large r is F[r) ~ C jr^ with a constant C. 
In the chiral bag, the pion carries a fractional baryon number 
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which is characterized by a topological property. By substituting the hedgehog solution, 
the fractional baryon number is given by 
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where F is a chiral angle at the bag surface |^. 

In the ud quark sector, we consider a basis set for the quark wave function by using a 
hedgehog solution h. We construct the quark wave function ip^'^'> for natural {k = +) and 
unnatural {k = — ) parity assignment for the grand spin K, respectively; 
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Jk{p^) is a spherical Bessel function with the ud quark momentum p. Here the basis 
states are given by 

|0) = YKMi9A)Xo, (18) 

|1) = Yl {KM - filfi\KM)YKM^^{e,<l))xl 

^=-1,0,1 

|2) = Yl (K +1M - filfi\KM)YK+iM-,i0A)xl, 

M=-1,0,1 

|3) = Y {K -IM - ^lfi\KM)YK-iM-,{0A)xl, 

^=-1,0,1 

where Yxm{0,4') is a spherical harmonics with spherical coordinate {0,(1)). Xu are eigen- 
states of the sum of spin and isospin G = S + I = a/ 2 + r/2, 

x^ = ^(IT)M)-|i)l«)), (19) 

xi = IT)H, 

xl = ^{\W) + \i)\u)), 

X\ = \i)\d). 

The sign of naturalness and unnaturalness corresponds to the parity P = {—)^^'^. The 
coefficients Oj and bi (i = 0, ■ ■ ■ , 3) are determined by satisfying the equation of motion of 
the ud quark. 

The eigenenergies of the hedgehog ud quarks are given by the boundary condition at 
the bag surface r = R, 

in ■ 77/'('') = -e^^-"-^(^)T5^('') , (20) 

which represents the continuity of the vector current flux. Then, we obtain the eigenvalue 
equation which is explicitly shown by 
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with the quark energy E = Jp^ + mj^. The equation for i^ = is obtained by setting 
j-i{pR) = 0. We note that there is a symmetry between energy levels in the positive and 
negative energy sides, 

Ej,p{F) = -Ek-p{-F). (22) 



This relation stemmed from the invariance of the lagrangian (^ under the transformation 
U ^U* oi F ^ -F. 

In the chiral bag model, the vr meson cloud at the bag surface induces the vacuum 
polarization in the quark sector, and consequently the change of the baryon number 
carried by the vacuum quarks. In the conventional chiral bag model with massless quark, 
it has been known that the fractional baryon number of the vacuum quarks is canceled by 
that of the pion cloud [^niE^- We have checked that this is also the case for the massive 
quarks. Namely, the fractional baryon number of the vacuum quark, 

B,{m,F) = -^lhn5^sgn(K)e-*l^"l (23) 

= -(f- -sin 2F 
vr V 2 

is independent of the quark mass in an analytical method by using the Debye expansion 
developed in [Hj. We also have checked that this result is supported by the numerical 
computation by the Strutinsky method. The detail procedure of the Debey expansion 
and the Strutinsky method for the massive quark will be discussed in the forthcoming 
paper. 

Since the classical hedgehog solution is not an eigenstate of spin nor isospin, the states 
obtained above are not physical states. In order to obtain physical states with proper 
quantum numbers, we need to quantize the classical hedgehog configuration through ro- 
tation in flavor SU{3) space. In this paper, however, we consider only the classical 
hedgehog solutions. As we will see, the strength of the pion field is not very strong, where 
it is known that the effect of the rotation is not important for the mass of a quark droplet. 

2.3 Energy in the NJL chiral bag 

In this subsection, we discuss the energy of a quark droplet including the energies of 
quark vacuum, valence quarks and the pion cloud. In order to estimate the quark vacuum 
energy, we consider the following two components. One is the energy caused by the NJL 
interaction which is expressed by an effective bag constant B and the other is the chiral 
Casimir energy. 

First, we evaluate an effective bag constant as the difference of the vacuum energies 
between the inside of the bag and the physical vacuum outside of the bag. The effective 
bag constant is obtained by considering the sea quarks, which are interacting through the 
NJL interaction in the bag. This is an empty bag in which sea quarks within the cutoff 
momentum contribute. It is then measured from the vacuum energy of the bulk vacuum 
in the NJL model. In the previous investigations J221I2S1; it was shown that the effective 
bag constant in this procedure agreed with the phenomenological value which is used in 
the MIT bag model. The effective bag constant is given by 
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where V = (47r/3)i?^ is the volume of the bag, N^ the number of the degree of freedom 
of color. Here K^ refers to the eigenstate in the vacuum part in the bag for q = u, d 
and s, respectively. The vacuum quark contribution is cut off by a smooth regularization 
function of a Lorentzian type ,23^ with a momentum cutoff A 

g{p/A) = -— — (25) 

^' ' l + (p/A)'^ ^ ' 

for smearing of the discrete energy levels of the quarks in the bag. The smoothness 
parameter a in the regularization (j^K|) is fixed to reproduce the nucleon mass E^+a = 1-1 
GeV. We have obtained a = 22.58. The last term is the energy in the bulk vacuum eg, in 
which chiral symmetry is maximally broken. We mention that the effective bag constant 
is given at the zero chiral angle. The effect of the pion cloud in the energy of the quark 
droplet is accounted by the chiral Casimir energy as discussed below. 

The chiral Casimir energy has been studied extensively in 80's and the result is well 
known. It is defined by 

Ec{m,R,F) = -i lim ^sign(E„)^„e-*l^"l. (26) 

Here, the sum is taken over all the states with positive and negative energies. In the 
massless case, it has been known that there is a logarithmic divergence in ()2f)|l |^ |33] . 
In order to remove this divergence, we impose a condition that the second derivative 
of the chiral Casimir energy with respect to F vanishes at F = 0. We found that the 
logarithmically divergent term is proportional to sin^ F for a finite quark mass by using 
the Debye expansion. This analytical technique has been developed in the chiral bag 
for massless quarks jUj. Subtracting the divergent term, we find the subtracted finite 
Casimir energy. 
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The reference point of the Casimir energy is set at F = 0. In order to calculate the sum 
numerically, the regularization, such as the exponential type, Gaussian type, heat kernel 
type |inj and so on, have been used. In our discussion, we use the Strutinsky smearing 
method |1I1I12J. This method has an advantage that the maximum energy in the sum of 
()2(i|l to achieve good convergence can be taken relatively as small as ERmax ~ 40, while 
at least ER^ax ~ 100 is necessary for other regulators. We note that ER^ax ~ 40 is 
needed for the case of massive quarks, while even a smaller value ER^ax ~ 10 is sufficient 
for massless quarks. Also, we note that the chiral Casimir energy always takes a positive 
value for any quark mass and chiral angle. 

Consequently, the total energy of the NJL chiral bag is given by a summation of the 
valence quark energy, the effective bag constant, the chiral Casimir energy and the pion 
cloud energy, 

Etotim, R, F) = E,ai{m, R, F) + B^ffim, R)V + El^^'im, R, F) + E^{R, F), (28) 



where m is the dynamical quark mass of either the ud or s quark. The first term indicates 
a contribution from the valence quarks. We assume iV„ = N^ in the discussion below, and 
further ignore the Coulomb energy, since we consider only small A system. 

The input parameters in the calculation of the energy of a quark droplet are the baryon 
number A and the strangeness S. For a given baryon number A and a strangeness S, 
Nu + Nd = NcA + S and A^^ = —S. The energy is then obtained by taking a variation 
with respect to the dynamical quark mass m^ = m^ and m^, the chiral angle F at the 
bag surface, and the bag radius R. 

3 Numerical Result 

3.1 Quark energy level in the NJL chiral bag 

The energy levels of ud quarks in the chiral bag are affected by the chiral angle at the bag 
surface F. By solving the boundary condition Eq. ()21|). we obtain the energy levels which 
are shown as functions of F in Fig. QJ The four figures correspond to the quark mass 
mR = 0, 1, 2 and 3, respectively. The K-grand spins are denoted for i^ = 0, 1, 2, 3 and 
4. The solid and dashed lines are for the positive and negative parity states. Horizontal 
dot-dashed lines are for \E\ = m. 

It is one of the particular points in the chiral bag model that the first 0^ states dives 
into the vacuum as the absolute value of the chiral angle F increases. The critical chiral 
angle for E = in the 0+ state is F/n = -0.5, -0.70, -0.82 and -0.87 for mR = 0, 1, 
2 and 3, respectively. As the quark mass increases, a larger chiral angle is necessary for 
pulling the valence quark down into the vacuum part. For O"*" state, we obtain E = m at 
F/tt = —0.5 for any quark mass. The states in \E\ < m are obtained by substituting the 
momentum with an imaginary value ip in Eq. (jUJ). As the quark mass increases, more 
states appear in the energy region of |-E| < -m, or \F\ ~ vr. For example, we have only the 
lowest 0^ states in this region for mR = 1, while 0^ and 1^ states are also included in 
|£^| < m for m,R = 2. For mR = 3, 2^ states are also included in \E\ < m. 

In the massless fermion, there is a symmetry of the energy level for F —>■ F + tt under 
the change in the parity. However, the massive fermion does not have this symmetry. In 
order to see the origin of this asymmetry, we write the equation of motion for the quark 
as 



mh - le*"-"^^^(^)(5(r - R) 



h = 0, (29) 



where the boundary condition is included as the delta function term. The sign of the 
5-function term changes from negative to positive as the chiral angle changes from F = 
to — TT. Therefore, the sum of the mass term and the scalar term proportional to cos F in 
the function changes from m^ + 6{r — R)/2 to m,h — S{r — R)/2 for F = to — vr. This 
is the origin of the asymmetry of the quark levels. These quantities can be interpreted as 
an effective mass, or the scalar potential, of the quark in the chiral bag. 
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Figure 1: The energies E of various quark levels in a chiral bag as functions of the chiral angle 
-F. (a) is for niR = 0, (b) mR = 1, (c) niR = 2 and (d) niR = 3, where m is the ud quark mass 
and R the bag radius. The horizontal dot-dashed lines indicate \E\ = m. The i^-grand spins 
are denoted for i^ = 0, 1, 2, 3 and 4. The it parity is distinguished by the solid [P = +1) and 
dashed (P = — 1) lines. 
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3.2 Multi-quark droplets 

Now, we discuss the stability of quark droplets. We are interested in the quark droplets 
with small baryon numbers up to A ~ 5, in which the surface pion is expected to have a 
non-negligible contribution. To start with, we discuss the procedure of energy variation 
of a quark droplet. We have several dynamical variables in the energy (j^H|) of a quark 
droplet with baryon number A and strangeness S; the dynamical quark mass (m^ = m^ 
and rris), the chiral angle F at the bag surface and the bag radius R. As shown in our 
previous work in the NJL model with the MIT bag [2Hj, variation of the total energy with 
respect to the dynamical quark mass leads to the restoration of chiral symmetry for -R < 6 
fm for ud quark sector {niud = 0), and i? < 3 f m for the s quark sector (m^ = mos). This 
is also the case in the NJL model with the chiral bag. 

A = l 

First, we discuss a strangelet of the baryon number A = 1, namely the nucleon. As 
mentioned above, since the chiral condensate (qq) vanishes, the results of the present 
model differs from the conventional bag model only by the bag constant term. The stable 
solution for a given bag radius R is then obtained by taking a variation of the total energy 
with respect to F. We show the numerical result in Fig. |2fa), where the total energy as 
well as various contributions, such as the pion cloud energy, the quark energy (the sum 
of the valence and Casimir energies), the Casimir energy and the volume energy of the 
effective bag constant, are shown as functions of R. The total energy has two local minima 
at i? = and 0.750 fm. The first local minimum corresponds to the pure Skyrmion, in 
which the quark contribution vanishes. The second minimum is characteristic in the 
present model, where both the pion and quarks contribute to the physical quantities. 
The chiral angle is shown as a function of the bag radius by the solid line in Fig. El At 
R = 0.750 fm, we obtain F/n = -0.327. 

The appearance of the second local minimum is due to the R dependence of the 
effective bag constant -Be// 1211 ■ Here, we note that the minimum value of the volume 
energy B^ffV depends on the smoothness parameter a in the regularization, while the 
corresponding bag radius does little. We have used that dependence in order to reproduce 
the mass of the nucleon. At this point, the resulting volume energy takes either positive 
or negative value. (In the present model, it takes a negative value.) However, its absolute 
value is not large, and the most part of the mass of the nucleon is given by the quark and 
the pion energies; E^ = 0.437 GeV and E^ = 0.715 GeV. 

A = 2 

Next, let us discuss the case of baryon number A = 2. We consider several configurations 
for ud quarks for a fixed strangeness S. The ud quarks can occupy various levels associated 
with the chiral angle as shown in Fig^a). The first three quarks occupy the 0+ state, 
and form a closed shell. For the remaining three quarks, the 1^ state is available for 
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F/tt > —0.317, while the 1 state is also used for —1 < F/n < —0.317. Therefore, we 
shall compare the following two configurations. 

For ^ = 0, we consider a) (0+)3(l+)3 and b) (0+)^ (1-)^. Here (K^)^ denotes a 
configuration how A^ quarks accommodate the state K^. The state K^ can contain 
Nc {2K+1) ud quarks maximally. By comparing the minimum energy of the configurations 

a) and b), we find that the configuration a) is more stable. The reason that we consider 
two configurations a) and b) is that lower quark levels depend on the chiral angle. We 
also note that the parity affects for a) and b). Here, we simply find a configuration which 
has the minimum energy. 

In Fig. Efb), we show the total energy as well as various components such as the 
quark energy (the valence quark plus the chiral Casimir energy), the Casimir energy, the 
pion cloud and the volume energy for ud quark configuration (0"'")^(1^)^ . We find a 
local minimum with the energy E = 3.042 GeV at the bag radius R = 0.750 fm, where 
the total energy is dominated by the valence quarks. They are £^3/1,0+ = 1-160 GeV and 
-E'3h,i+ = 1-770 GeV for the 0+ and 1+ states occupied by three quarks, respectively, 
while the pion cloud energy is E^r = 0.148 GeV which is smaller than that in the nucleon 
{A = 1). The small pion contribution is due to the small chiral angle in absolute value. 
As shown in Fig. El we find F/n = —0.165 at the stable point, which is smaller than that 
of the nucleon F/n = —0.327. 

For a finite strangeness, we investigate the configurations a) (0+)^ (i+)3-|5| ^^^^^ ^-^ 
(O"*")^ (l^)^~l'^l for S = —1 and —2, respectively. Again, the configuration a) is found 
to be more stable for each strangeness. We summarize the numerical results in Table ^ 
where we show the ud quark configurations, the energy per baryon number E/A, the bag 
radius R, the chiral angle F and the pion energy E^^ for each strangeness. The energy per 
baryon number is E/A = 1.521, 1.527 and 1.526 GeV for S = 0, —1 and —2, respectively. 
Therefore, the most stable state is the non-strange state of 5* = 0, when the bag radius is 
R = 0.750 fm, and the chiral angle F/vr = —0.165. 

To see the reason that the non-strange state is more stable, let us look at the energy 
levels of Fig. 1(a). At the chiral angle F/tt = —0.165 where R = 0.750 fm, among the six 
ud quarks, the three ud quarks occupy the O'^ state, and the other three occupy the 1+ 
state with Eh^i+ = 2.24/ R. The latter should be compared with the energy of the strange 
quark of 1/2+ state, £'5,1/2+ = 2.30/R. Hence s quark state is above the ud quark 1+ 
state, and therefore, the strangeness cannot be included in the ground state. 

From the above discussions, we have seen that the energy per baryon number E/A = 
1.521 GeV of the quark droplet is larger than the nucleon mass -Eat+a = 1.1 GeV. There- 
fore, our result indicates that the quark droplet with A = 2 is not stable against the decay 
into two nucleons. 

^ = 3 and 4 

For y4 = 3, we investigate the following ud quark configurations, a) (0"'')'^ (i+)6-|S'| ^^^^ 

b) (0+)^ (i-)6-|S'| fQj, ^ _ g^ _^^ _2 g^]2(i _3^ and found that the configuration a) is 
more stable than b) for each strangeness. We show the numerical results in Table |2l 
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Figure 2: Various energy contributions to a quark droplet of baryon number A = 1 (nucleon) 
and A = 2. 
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Figure 3: The chiral angle F at the bag surface as functions of the bag radius R for non-strange 
quark droplets of the baryon number A = 1, 2, 3 and 4. 
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We obtained the energy per baryon number E/A = 1.588, 1.605, 1.618 and 1.630 GeV 
for S = 0, —1, —2 and —3, respectively. Therefore, the most stable state is the non- 
strangeness state with S = 0. There, the bag radius is R = 0.754 fm, and the chiral 
angle is F/n = —0.0580. The chiral angle is very small as compared with that of the 
nucleon. Correspondingly, the pion cloud energy is E.„- = 0.02267 GeV, which is very 
small as compared with that of the nucleon. 

In order to understand the reason that the strange quark is not contained in the 
ground state, we turn to the ud quark energy level in Fig. ^a). At the chiral angle 
F/n = —0.0580, the three ud quarks occupy the lowest 0"*" state. The other six ud quarks 
occupy the 1+ state with eigenenergy Ei+ = 2.24//?, which is lower than the lowest s 
quark energy i^^ 1/2+ = 2.30/ R. The 1~ energy is also lower than the s quark energy, and 
therefore, the s quark does not appear in the ground state of the quark droplet for A = 3. 

For the baryon number A = A, we considered the configurations a) (0^)^ (^i+)9-l'S'l g^Yid 
b) (0^)^ (1^)^^'"^' for the strangeness 5 = 0, —1, —2, —3 and —4. We have found again 
that the configuration a) is more stable for each strangeness. We show the numerical 
result in Table El where we see that the energies per baryon number are E/A = 1.597, 
1.616, 1.630, 1.644 and 1.731 GeV for S = 0, -1, -2, -3 and -4, respectively. Hence, 
the non-strange state of 5* = is the most stable one, when the bag radius is R = 0.757 
fm and the chiral angle is F/n = —0.144 x 10~^. Here, the chiral angle is negligibly small, 
and so is the pion cloud energy. In this way, it is justified to use the MIT bag model 
instead of the chiral bag model for A = A and larger fl^ ESj (see also the discussion 
below) . 

A = 5 

Lastly, we consider the case of baryon number A = 5. We consider the ud quark configu- 
rations a) (0+)3 (l+)9 (l-)3-l^l and b) (0+)^ (l-)^ (l+)3-|5l and c) (0+)^ (l-)^ (O-)^-!^! for 
S = 0, -1, -2, and a) (0+)^ (l+)9-l^l and b) (0+)^ (l-)9-|5| fo^ S = -3, -4, -5. Among 
them, the configuration a) is the most stable one for each strangeness. In Table IH we 
show the numerical results. The energy of the quark droplet is E/A = 1.762, 1.725, 1.683, 
1.638, 1.710 and 2.226 GeV for S* = 0, ■ ■ ■ , —5, respectively. Therefore, the most stable 
state is S* = — 3 of finite strangeness at the bag radius R = 0.761 fm, where the chiral 
angle and the pion cloud energy are almost zero. 

We can clarify the reason why the finite strangeness is included in the ground state of 
the quark droplet with A = 5. Let us turn again to the ud quark energy level in Fig. ^a). 
At the chiral angle F ~ 0, the three quarks occupy the O"*" state with eigenenergy Eh^Q+ ~ 
2.04/i?, and the nine quarks form a closed shell in the l"*" state with Efii+ ~ 2.04/i?. The 
remaining three quarks may occupy the 1~ state with E^^i- ~ 3.20/-R, which is larger 
than the strange quark energy -Es,i/2+ = 2.31/i?. Therefore, the three quarks occupy the 
1/2"'" state in the s quark sector. Concerning the chiral angle F, it is exactly the same 
for the case oi A = 4, since the ud quark configurations are the same. Therefore, it is a 
good approximation to neglect the pion cloud for multi-baryon system. 

For the baryon number A > 5, the results of the NJL chiral bag model becomes 
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similar to the N JL model with the MIT bag jzzl 123] . As discussed in our previous work 
|23j , the Multiple Reflection Expansion (MRE) reproduces the energy per baryon number 
obtained by solving the discrete quark levels for the baryon number A > 10. There, it 
was concluded that the strangelets with baryon number 10 < A < 10^ were not stable as 
compared with the normal nuclei P2j . In this paper, we have verified again that the N JL 
chiral bag model does not support a stable strangelet as the ground state of the finite 
quark system. This conclusion as derived in our previous works is qualitatively different 
from the analysis by the MIT bag model [3 El E] • 



s 


/j^P\#{ud quarks) 


E {E/A) [GeV] 


R [fm] 


F/tc 


E^ [GeV] 





(o+f(i-f 


3.042 (1.521) 


0.750 


-0.165 


0.1480 


-1 


mHi-y 


3.054 (1.527) 


0.750 


-0.197 


0.2117 


-2 


{o^ni-y 


3.052 (1.526) 


0.749 


-0.230 


0.2856 



Table 1: The ud quark configurations, the total energy E, the bag radius R, the chiral angle F 
and the pion cloud energy £J^ for quark droplets A = 2 with several strangeness S. The values 
in the parentheses in the third row are the energies per baryon number E/A. 



s 


/ Ty-P\#(ud quarks) 


E {E/A) [GeV] 


R [fm] 


F/'K 


E^ [GeV] 





(0+)^!+)^ 


4.764 (1.588) 


0.754 


-0.0580 


0.02267 


-1 


(o+)ni+)5 


4.815 (1.605) 


0.754 


-0.0811 


0.04429 


-2 


{o^ni^r 


4.854 (1.618) 


0.754 


-0.105 


0.07527 


-3 


{o^ni^r 


4.890 (1.630) 


0.753 


-0.131 


0.1169 



Table 2: Same as Table 1 for A = 3. 



s 


/j^P\#{ud quarks) 


E {E/A) [GeV] 


R [fm] 


F/tt 


E^ [GeV] 





(0+)^!+)^ 


6.388 (1.597) 


0.757 


-0.144 X 10"^ 


0.1675 X 10"^ 


-1 


(0+)^!+)^ 


6.464 (1.616) 


0.750 


-0.0151 


0.001828 


-2 


(0+)^!+)^ 


6.520 (1.630) 


0.757 


-0.0308 


0.007702 


-3 


(0+)3(l+)6 


6.576 (1.644) 


0.757 


-0.0481 


0.01876 


-4 


(o+)ni+)^ 


6.924 (1.731) 


0.758 


-0.0667 


0.03612 



Table 3: Same as Table 1 for A = 4. 



4 Conclusion 

In this paper, we discussed the stability of strangelets by considering dynamical chiral 
symmetry breaking of QCD. We investigated the effects of the dynamical generation of 
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s 


/T^p\^{ud quarks) 


E {E/A) [GeV] 


R [fill] 


F/tt 


K [GeV] 





m^{i^f{i-f 


8.630 (1.762) 


0.762 


-0.119 


0.09780 


-1 


r)niT(i')^ 


8.625 (1.725) 


0.762 


-0.0702 


0.03348 


-2 


m^{i+f{i~Y 


8.415 (1.683) 


0.761 


-0.0200 


0.03820 


-3 


(o+)ni+r 


8.190 (1.638) 


0.761 


-0.854 X 10-^ 


0.6923 X 10-^ 


-4 


m^{i^f 


8.550 (1.710) 


0.761 


-0.0127 


0.001536 


-5 


m^{i^y 


11.13 (2.226) 


0.762 


-0.0304 


0.007548 



Table 4: Same as Table 1 for A = 5. 

quark masses in a finite volume by introducing the NJL model inside the chiral bag. This 
is a chirally symmetric two phase model, in which quarks inside the bag are interacting 
through the point-like interaction of the NJL model, and mesons exist outside the bag as 
meson cloud. 

We adopted the hedgehog configuration to calculate the energies of the systems of 
baryon number A < 5. The dynamical quark masses of the m, d and s quarks were then 
given in the mean field approximation of the NJL interaction. It turned out that for 
strangelets with small baryon number, chiral symmetry is restored inside the bag, and 
hence the dynamical mass of the quarks vanishes. However, for general purposes, we have 
investigated various chiral Casimir effects with finite quark mass. By using analytical and 
numerical methods, it was verified that the fractional baryon number in the ud quark 
sector is not affected by the finite quark mass. Therefore, the total baryon number is 
conserved exactly, as in the case of massless quarks. The Casimir energy was also well 
defined by subtracting the divergent term for the massive quark. 

The total energy of a quark droplet was expressed as a sum of the energy of the 
valence quarks, the Casimir energy, the volume energy of the effective bag constant and 
the pion cloud energy. The energy of the system was measured from the reference point 
of the chirally broken vacuum outside the bag. The energy of the quark droplet with 
baryon number A was determined by taking a variation with respect to (in addition to 
the dynamical quark masses) the bag radius, the chiral angle and the strangeness. As 
a result, we could reproduce the nucleon mass E = 1.1 GeV and the nucleon radius 
R = 0.750 fm. For the multi-baryon states, it was shown that the quark droplets with 
baryon number A = 2, 3 and 4 did not contain the strangeness in the ground state. 
For the baryon number A > 5, the quark droplets contained the finite strangeness. We 
obtained the strangeness S = —3 for A = 5. For complication of analysis of discrete 
states, we did not investigated explicitly systems oi A > 5, but we expect that for larger 
A strangelets of finite strangeness will become the ground states. 

We introduced the Lorentzian regularization for the quark vacuum energy associated 
with the NJL interaction. This quantity was expressed by an effective bag constant. We 
have also investigated the Gaussian type regularization instead of the Lorentzian. We 
found that the Gaussian regularization gave the effective bag constant which was larger 
than the value used in the conventional MIT bag model, which is consistent with our 
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previous result PHJ. Accordingly, the mass of the nucleon became too large and could not 
reproduce the observed value. 

In the present study, we have found that the pion contribution was suppressed as the 
baryon number increased. The pion cloud energy in the quark droplet oi A = 5 (and 
likely for larger A) was almost negligible. In this way, the use the MIT bag model for 
the quark droplets with baryon number A > 5 |221 llBj can be justified. The role of the 
pion cloud was previously investigated for non-strange multi-baryon systems with A = 2, 
3 and 4 in chiral bag model ^7\ EH IHH HO] ■ Our result is consistent with theirs. 

The energies per baryon of these quark droplets are larger than the mass of the nucleon. 
Therefore, the strangelets are not stable against decays into multi-baryons. This result is 
qualitatively different from the previous studies in the MIT bag model. 

There are several issues to be further discussed. In this paper, we have assumed 
the hedgehog ansatz at the classical level. For realistic comparison with experimental 
data, it is necessary to perform spin-isospin projection. It is also interesting to discuss 
the Nambu-Goldstone mode in the chiral symmetry broken phase in the bag for large A 
system. There, the Nambu-Goldstone mode will directly couple to the quarks inside the 
bag, not only at the bag surface. The vacuum structure in the bag may be affected by 
this volume type coupling. We are now currently working on this problem. The explicit 
U{1)a breaking is also an important topic. There, the 't Hooft term will be introduced in 
addition to (or instead of) the NJL interaction. The effect of the qq correlation is also an 
interesting subject ^HlElllHl- As the baryon number increases, the number of the valence 
quarks also increases. There, it would be possible to cause the color superconductivity. 
It is an interesting question whether the qq correlation in the bag is affected by the pion 
cloud outside the bag. 
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